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We examine the coexistence of the vortex state and the Fulde-Ferrell-Larkin-Ovchinnikov 
(FFLO) state in quasi-two-dimensional type-II superconductors and the crossover from the coexis- 
tence state to the pure FFLO state when the Maki parameter a increases. The pure FFLO state, 
characterized by finite center-of-mass momenta q 7^ of Cooper pairs occurs in the two-dimensional 
limit, when the magnetic field is parallel to the conductive plane. The vectors q are determined from 
the Fermi-surface structure and pairing anisotropy, and become finite below a temperature T* . In 
quasi-two-dimensions, because of the orbital pair-breaking effect, the coexistence state characterized 
by (n, Qii) occurs, where n and denote the Landau level index of the vortex state and the wave 
number of the additional FFLO modulation along the magnetic field. We obtain the a dependence 
of the upper critical field by numerical calculations. The upper critical field exhibits a cascade curve 
in the H-T phase diagram. It is analytically shown that n diverges in the two-dimensional limit 
Q — > 00 below T* . In this limit, the upper critical field equation of the coexistence state is reduced 
to that of the FFLO state. A relation between n of the coexistence state and q± of the pure FFLO 
state is obtained, where q± denotes the component of q perpendicular to the magnetic field. It is 
found that the pure FFLO state is nothing but the vortex state with infinitely large n as is known 
in two-dimensional superconductors in a tilted magnetic field. The vortex state with large n can be 
regarded as the FFLO state with non-zero q± in three dimensions. 



PACS numbers: 74.81.-g 74.25.-q, 74.25. Op, 
I. INTRODUCTION 

In type-II superconductors, an applied magnetic field 
destroys the superconductivity by two kinds of pair- 
breaking effects: the orbital magnetic and Pauli param- 
agnetic pair-breaking effects ,1] . We can define the pure 
orbital limit i?c20 and the Pauli paramagnetic limit Hp, 
which are theoretically obtained by taking into account 
only the orbital effect and paramagnetic effect, respec- 
tively. The Maki parameter a = y/2Hc2t)/ Hp expresses 
the strength ratio of the two pair-breaking effects. 

In conventional metal superconductors, the orbital 
pair-breaking effect dominates the system (i?c20 ^^p) 
because of the large Fermi velocity. Partial destruction 
of the superconductivity due to the orbital effect creates 
vortexes, which form a lattice below the upper critical 
field, and causes the order parameter to become nonuni- 
form. 

In contrast, in purely Pauli limited superconductors, 
another type of nonuniform superconductivity has been 
proposed by Fulde and Ferrell [3| and Larkin and Ovchin- 
nikov 3. In magnetic fields, the Fermi surfaces of the 
up- and down-spin electrons are displaced due to the Zee- 
man energy. If the up- and down-spin electrons on the 
displaced Fermi surfaces form Cooper pairs, they should 
have a finite center-of-mass momentum. The supercon- 
ducting state of such Cooper pairs is called the Fulde- 
Ferrell-Larkin-Ovchinnikov (FFLO) state. It is easily 
verified that the finite center-of-mass momentum results 
in spatial modulations and nodes of the order parame- 
ter in real space. As a result, spin polarization energy 



is gained by the depaired electrons near the nodes, while 
the condensation energy is lost. 

Therefore, a necessary condition for the occurrence of 
the FFLO state is that the superconductivity survives in 
high fields such that ^^H ^ Aq, for which the loss of 
condensation energy due to modulation of the order pa- 
rameter is compensated by a gain in polarization energy. 
Here, /ie and Ag denote the electron magnetic moment 
and the zero field energy gap of the superconductivity. 
This condition can be expressed as Hc2 ~ Hp, where 
Hc2 denotes the upper critical field, because Hp ^ Aq. 
Therefore, the orbital pair-breaking effect needs to be 
very weak for the FFLO state to occur. 

This is one of the reasons why the FFLO state has not 
been observed in conventional metal superconductors. 
Gruenberg and Gunther found that the FFLO state oc- 
curs only when the Maki parameter a is large (a ^ 1.8) in 
isotopic superconductors Such a large Maki param- 
eter is usually difficult to realize in alloy type II super- 
conductors, though it is achievable in some exotic super- 
conductors, such as organic superconductors Q, heavy 
fermion superconductors Q, and oxide superconductors, 
because of their narrow electron bands, large effective 
masses, and quasi-two-dimensionality. 

Because of the orbital pair-breaking effect, and the for- 
mation of the vortex lattice state, the dependence of the 
order parameter on the spatial coordinates perpendicu- 
lar to the magnetic field is described by the superposi- 
tion of the Abrikosov functions. Therefore, there may 
be additional modulations due to the FFLO state only 
in the direction parallel to the magnetic field, as Gruen- 
berg and Gunther have proposed We examined the 
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coexistence of the vortex states with higher Landau level 
indexes n and the FFLO state in d-wave superconduc- 
tors in our previous paper [3], and obtained phase dia- 
grams for some cases. Recently, in a model of the quasi- 
two-dimensional heavy fermion superconductor CeCoIns, 
Adachi and Ikeda have shown that the first-order phase- 
transition to a coexistence state occurs taking into ac- 
count the higher Landau level indexes @. 

Bulaevskii examined film superconductors at T = 
in tilted magnetic fields, and obtained cascade transi- 
tions between the states with different n's when the di- 
rection of the magnetic field changes. It was shown 
that the upper critical field tends to approach that of 
the FFLO state when the magnetic field approaches the 
parallel direction. Buzdin and Brison also obtained cas- 
cade transitions at finite temperatures in two- and three- 
dimensional isotropic superconductors @. We extended 
Bulaevskii's theory to finite temperatures in s- and d- 
wave superconductors [13] ■ We also reproduced the cas- 
cade transitions solving the gap equation, and clarified 
the behavior in the limit of a parallel field. We analyti- 
cally demonstrated that the Landau level index n of the 
vortex state diverges in the temperatures region where 
the FFLO state occurs in the limit. As a result, the 
envelope of the cascade transition lines approaches the 
FFLO critical field when the field orientation approaches 
the parallel direction. The coexistence state is continu- 
ously reduced to the pure FFLO state. We obtained a 
relation between n and the FFLO vector q, which con- 
nects the vortex states and the FFLO state in the limit. 
The present study extends our previous theory to the 
coexistence state in three dimensions. 

Near critical fields, the order parameter of the pure 
FFLO state can be expressed by a linear combination 
of exponential functions exp[iq„i • r], where are the 
degenerate FFLO wave vectors with optimum values de- 
termined from the structures of the Fermi surfaces and 
the pairing interactions, and the temperature. The free 
energies of such states were compared in a three dimen- 
sional isotropic system by Larkin and Ovchinnikov 
Matsuo et al. Illll. B owers and Rajagopal [l^, and Mora 
and Combescot [IJ], and in two-dimensional systems by 
the author [T^ and Mora and Combescot [3]. Larkin 
and Ovchinnikov [sj showed that the state expressed 
by A(r) ~ cos(q • r) has a lower free energy than the 
state A(r) ^ exp[iq • r] proposed by Fulde and Ferrell. 
The author has shown that the square, triangular, and 
hexagonal states have lower free energies than the state 
A(r) ~ cos(q • r) at low temperatures in two dimen- 
sions jljj. Mora and Combescot have shown that the 
states described many cosine functions have lower free 
energies than the state A(r) ~ cos(q • r) in three dimen- 
sions when the first-order phase transition is taken into 
account [T^ . The possibility of the first-order transition 
was addressed also by Larkin and Ovchinnikov Q . 

When the magnetic field is not oriented in the optimum 
direction of q of the pure FFLO state, or when more 
than two qm's contribute to the linear combination for 



the pure FFLO state, it may appear that the coexistence 
state is not reduced to the pure FFLO state in the limit 
a oo, because in the coexistence state q can have a 
nonzero component only in the direction of the magnetic 
field. In actuality, however, the Landau level index n 
of the coexistence state diverges and the components of 
q perpendicular to the magnetic field are realized ^16j. 
This behavior is analogous to that in two-dimensional 
systems in a tilted magnetic field. In the present paper, 
we demonstrate this behavior in quasi-two-dimensional 
systems by an analytical proof and concrete numerical 
calculations. 

For simplicity, we adopt an effective mass model 
assuming that the first-order transition is suppressed, 
to demonstrate the continuity between the coexistence 
states and the pure FFLO state. In actuality, in the ef- 
fective mass model with very large a, the first-order tran- 
sition would occur at a slightly higher field than that of 
the second-order transition [13]. In this case, the second- 
order transition curve below the first-order transition 
curve is regarded as that of the metastable transition, 
which can be realized when the system is supercooled. 

Recently, we have obtained a result which also sug- 
gests that the Landau level index n increases with a in 
an anisotropic Ginzburg-Landau model near the tricrit- 
ical point [ITj ]. Comparing the phase diagrams with and 
without the orbital effect, we have found that the areas 
of the coexistence states with n > in the former phase 
diagrams correspond to the areas of the pure FFLO state 
with q± m the latter phase diagrams. 

In section 2, we present our formulation. In sections 
3 and 4, we examine the pure FFLO state and the co- 
existence state, respectively. In section 4, we show the 
crossover from the pure FFLO state to the coexistence 
state and the limit a ^ oo. In section 5, we show the 
numerical result for finite temperatures. In section 6, we 
summarize and discuss the results. 



II. FORMULATION 

We examine a model described by the Hamiltonian 

H = Ha + H,n + H', (1) 

with 

Hrr, = Y, [d'rah^l{r)Mr), (3) 



(4) 
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Here, we have defined the effective masses mi = rrix, 
m2 = my, 7713 = ™zj Zeeman field h = /ie|-ff|, and vec- 
tor potential A, where fie and H denote the magnitude 
of the electron magnetic moment and the magnetic field 
H = rotA. We consider pairing interactions of the form 



V{p,p') = gala{p)la{p'), 



(5) 



where the suffix a expresses a symmetry. 

In the effective mass model of Eq. ([2]) , it is convenient 
to define f = (xi, 0:27 is) by a scale transformation 

y/mJ^Xf, = V^Xf_, (6) 

with m — {rnxmym^y^^ . Then, Eq. ([2]) is written as 

i^o^E/^^-'/^iwiH^J;-^^.)V.M, (7) 

/icr ^ 



where we have defined 




(8) 



We also define p = (pi,p2,P3) with = (m/mp)^/^^^, 
so that r ■ p — f ■ p. Then, the Fermi surface in p 
space becomes spherically symmetric, and we can de- 
fine a constant Fermi momentum pp and Fermi velocity 
vp = vfp/\p\ with a constant magnitude vp = pp/'fh, 
for the scaled momentum p. 

Now, we derive the gap equation. The calculation is a 
straightforward extension of the previous studies 0, fiol . 
[I8, 19]. Near the second-order phase transition, the gap 
function has a form 



A(r,p) = A„(r)7„(p), 
and the gap equation is linearized as 
T 



(9) 



log 



T, 



ttT 



(0) 



)A„(r) 



dt 



1 



X 



sinh(7rrt) J At: 



(10) 



l-cos[t{/i--i;^-n}]jA„(r), 



where v'. 



F - vfP'/\p'\ and n = (Hi, 02,03) with 
d 



c 



(11) 



The upper critical field Hc2 is the highest \H\ among 
those which give a nontrivial solution of A(r). We note 
that Eq. is the same as that of a system with a spher- 
ically symmetric Fermi surface except that the argument 
p' ~ p'/Ip'I of la{p') is different from the integral vari- 
able p — p' /\p'\ — p' Ipf- Therefore, when 7q,(p') is con- 
stant, it is easily verified that the mass anisotropy does 



not affect the upper critical field equation except that the 
vector potential is scaled as described in Eq. ([5]), since p 
is only an integral variable. In contrast, for anisotropic 
superconductors, the mass anisotropy affects the upper 
critical field equation through the deformation of ja{p') 
when it is expressed in p space [l^ . 



III. 



THE PURE FFLO STATE 



In this section, we briefly review the case in which the 
orbital pair-breaking effect is negligible. In this case, we 
can set A = in Eqs. ([2]) and (fTTj) . Equation (fTO)) has a 
solution of the form 



A(r) cx exp[iq • f/H], 



(12) 



and is reduced to 



log( 



^(0)J 



An 



\lc.{p'\ 



dt 



sinh(7rTi) 



1 - cos[t{h ~ -v'p ■ q]] 



(13) 

For example, for s-wave pairing, ^s{p) = 1, there is 
infinite degeneracy with respect to the direction of q, al- 
though the magnitude |q| is uniquely determined so that 
the critical field is maximized. For non-s-wave pairing, 
both the direction and the magnitude of q are optimized. 
Depending on the symmetries of 7q(p) and the Fermi 
surface, and the temperature, there may be 2, 4, 8, 16 • • •- 
fold degeneracies with respect to the direction of q. We 
write the optimum q's as (jm with m = 1, 2, • • • , M. Be- 
low and near the upper critical field, the order parameter 
is expressed by a linear combination: 



A(r)^EA, 



(14) 



Among the states of this form, the physical state is that 
with the lowest free energy. Every degenerate (jm does 
not necessarily appear in the linear combination of the 
physical state. The most well-known form is that ex- 
pressed by the linear combination of e"^ '' and e^"' '', i.e., 
A(r) c>c cos(q • r). For the second-order transition and 
s-wave pairing, this state is the physical state in the ef- 
fective mass model. 

Here, we note that (jm are not necessarily parallel 
to the magnetic field, when the orbital effect is negli- 
gible. The number and directions of the optimum q,„ 
which contribute to the physical state also depend on 
the structures of the Fermi surface and the pairing inter- 
actions ^0|, [lH , and the temperature [13, [H, IH, [2I] . In 
the present effective mass model, since the electron dis- 
persion becomes isotropic in p space as seen in Eq. (|13p , 
the Fermi-surface anisotropy does not remove the infinite 
degeneracy of the optimum qm for s-wave pairing. 
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IV. 



COEXISTENCE STATE 



In this section, we take into account both the orbital 
and paramagnetic pair-breaking effects. After deriving 
the upper critical field equation, we take the limit of a 
weak orbital effect. 

For a magnetic field H = {0,0, H), we define A = 
{—Hy, 0, 0) with an appropriate gauge. The scale trans- 
formation described above gives = —Hy with H = 
(fh/ ^rnxrny)H . We define boson operators by 



(15) 



with 



2\e\H 



(16) 



where = y/c/2\e\H = ^J^^J^tkH, which is of the 
order of the BCS coherence length when H ~ iJc2- 
The operator v'p ■ 11, which appears in Eq. PH)) . can be 
rewritten as 



n 



sin6i'(e"^'i?-He~"^'fyt) 



-\fmp cos 



4 

az 



(17) 



where d' and (p' denote the polar coordinates when z-axis 
is the polar axis. 

Because Eq. ([TUl) can be regarded as an eigen equa- 
tion with the eigenvalue — log(T/Ti'^'') and eigenfunction 
AQ,(r), our problem is reduced to finding the eigenfunc- 
tions with the highest eigenvalue. The solutions can be 
written in the form 



A(r) A(x,y) exp[ig3z/n], 



(18) 



with djdz in Eq. (fT7|) replaced by iq^jh. The gap equa- 
tion (flUl) can be rewritten as 



log( 



(0) 



)0(5,y) 



dt 



1 



sinh(7rTt) ] in 



(19) 



1 



1 - cos[t{h- -vpq^cos0' - Q] 4>ix,y), 



where we have defined 

Up sin 0' 

c 



(e'^r) 



e-'-^'ryt) 



(20) 



It is convenient to expand the function (j){x,y) by the 



Abrikosov functions 4>ii\x,y) defined by 



(21) 



where n = 0, 1, 2, 3, • • • , called the Landau level indexes, 
k is an arbitrary wave number, and c/jq''^ is the solution 
of 



f]4>i\x,y) = 0, 



which is expressed as 

^^'=Hi,y) = Ce"=^exp 



{y - iikf 



(22) 



(23) 



with j/fc = k£,fj and a normalization constant C . The 
function (j)n\x,y) is expressed as 



4>^n\i,y) = (-l)"Ce"=^i/„ 



V2- 



y-yk 



X exp 



{y-yk? 



(24) 



in terms of the Hermitc polynomial. The operators fj and 
fj^ and the Abrikosov functions 0^'^' satisfy the relations 



f]^ <t'n\x,y) = + 1 (t}-^l^{i,y) 
f]4>n\x,y) = ^/n (l)-^\{x,y). 
If we expand the eigenfunctions as 

oo 

A(r) = ^ A„(/)W(i,y) exp[iq^z/;i], 

n=0 



(25) 



(26) 



the gap equation (jlOp can be written as a matrix 
equation for the eigenvector with the vector elements 
Ao, Ai, A2, - • •, as 



T 



log(^)A„ -^Z?„„,A, 



(27) 



with 



ttT I dt- , 

/o sinh(7rTi) 



X I ^bdp'f I didyc^i'\i,y) 



l-cos[t{h~ ^vfQz cos 9' - C)] (l)nHx,y). 



(28) 



When h = 0, for s-wave pairing, the solution with n = 
and Qz = gives the highest upper critical field. In 
general, Abrikosov functions with different ri's can be 
mixed. 

The magnetic field \H\ appears both in the Hamilto- 
nians Hq and Hm- The field \H\ which originates from 
A in Hq is responsible for the orbital effect, and appears 
in the gap equation as a dimensionless parameter a^l-ff | 
with the coefficient defined by 



m 2|e| 



( V 



(29) 
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In contrast, the field \H\ included in the Zeeman field 
h in is responsible for the Pauli paramagnetic pair- 
breaking eS'ect. For this field, it is convenient to define 
the dimensionless parameter fie\H\/{2TTTc^^). The rel- 
ative strength of the Pauli paramagnetic pair-breaking 
effect to that of the orbital effect is expressed by the ra- 
tio of their dimensionless parameters 



|Jf|/(2rfO 



Me 



(30) 



The parameter z,„ is proportional to the Maki parameter 
a. For example, for s-wave pairing, numerical calcula- 
tions give amHc2o ~ 1.0372 and /ieiJp/Ao w 0.707107, 
and hence a ~ 7.39 x Zm- 

If we define am and Zm for the isotropic system by 



Zr 



{2\e\/c){vF/2nTc 
Me/(27rTi°^am), 



(0)x2 



we obtain 



j mxmy \ 
\ ml J 



1/6 



1/6 



(31) 



(32) 



with the Maki parameter of the isotropic system a = 
V2Hc2o/Hp. For the magnetic field parallel to z-axis, 
when rrix S> rriz we obtain Zm ^ Zm, i.e., a ^ a, which 
means that the system is strongly Pauli paramagnetic 
limited. 



V. CROSSOVER FROM THE PURE FFLO 
STATE TO THE COEXISTENCE STATE 

Now, we consider a quasi-two-dimensional system such 
that rrix 3> my,mz. In this case, because we have 
^ £,H ^ $0 from Eq. we may omit ( ^ 

{e'^'fi+e-^"^\f)/^H in Eq. ^ for finite n and n'. There- 
fore, if we can truncate the summation over n in Eq. (pS)) . 
the upper critical field equation is reduced to 



log( 



T . 

da 



^ ttT 



dt 



1 



4tt 



sinh(7rrt) 
- cos[t(/i - 



-Wi^g^cosfii')] 



(33) 

in the limit — > oo, which coincides with the up- 
per critical field equation for the pure FFLO state ex- 
pressed by Eq. (fT2|) with q = {0,0, Qz). The magnitude 
\q\ = \qz\ should be optimized so that the upper crit- 
ical field is maximized. We write the optimum value 
as qq. For example, for s-wave pairing, it is known 
that qo « 1.2 x 2h/vF at T = 0. Thus, we only have 
two states with the highest upper critical field, which 
have the FFLO vectors q = (0, 0,±go). This contra- 
dicts the fact that there are more than two e'^'" ''/'' in 
most exactly two-dimensional systems as mentioned be- 
low Eq. (|12p . Furthermore, below the upper critical field. 



the free energy is minimized by the state expressed by 
the linear combination of more than two e'^™ ''/'' at low 
temperatures [H [H [11 . 

This contradiction is due to the assumption that in- 
finitely large n's are negligible in the two-dimensional 
limit. In the limit m^; 3> my,mz, we need to consider 
Abrikosov functions with infinitely large Landau level 
index n. For large n's, the states with n ± 1 can be 
approximated by the state with n, which means that 
(/)„±i « e^"^°(j)n in Eq. ^E^, where e^"^° are arbitrary 
phase factors. Therefore, we may write 



V 



(34) 



in the gap equation. This procedure is analogous to 
those in the theory of Bose condensation, and in two- 
dimensional type-II superconductors in a tilted magnetic 
field [10|. Using Eq. §^ we obtain 



c 



/nvp sine 

V2iH 



■ cos{(p' - ipo). 



(35) 



Thus, the gap equation Eq. (fTO|) can be rewritten for 

rrix S> my, mz as 



T 



/■°° , 1 f d^f, 

ttT / dt / — ^ 

Jo smhinTt) J An 

1 - cos[t{h - ^v'p ■ q)] (t){x,y). 



Hlcip't (36) 



where 



Q= {qx,qy,qz) = (g_icos<^o,q'±sin(po,9z) (37) 



with 



_ . 1,2 ^ „~2 _ ^ _ . (38) 



11- = \/<lx+<ly = 



If we write the optimum n for each fixed a as n{a), we 
obtain 



= lim 



' 2n{a)a 



(39) 



If 7^ in the limit a — > oo, n(a) must diverge like a, 
i.e., n{a) ^ a. Equation (|36p coincides with the upper 
critical field equation of the pure FFLO state with q. 

Equation (|39[) is an essential equation which connects 
the pure FFLO state with the optimum q in the two- 
dimensional limit and the coexistence state with the op- 
timum n and qz in quasi- two-dimensions. From Eqs. (j38p 
and p9|) , we can see that if there are more than two op- 
timum q's in the two-dimensional limit, there must be 
coexistence states with different n's with close upper crit- 
ical fields in quasi- two-dimensions where m^ ^ my,mz. 
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For example, a d-wave superconductor with 



7d„ 



(40) 



at low temperatures exhibits a degeneracy in the so- 
lutions of Eq. dM]) with q = {0,±qo,0), (0,0, ±qo) in 
the limit rrix — > oo. From Eqs. p7|) and (|39p . we 
find that when nix 3> my,mz, the coexistence states 

A(r) A„0l'=)(i,y)e'^-/'' with (n,g,) « (<z^el/2,0) 
and those with (n, q^) = (0, igo) have upper critical fields 
close to each other. 



VI. 



s-WAVE PAIRING 



In this section, we consider an s-wave superconductor 
as an example. Since 7s (p) = 1, the gap equation (fTUl) 
is exactly the same as that of the isotropic model except 
that the vector potential A is scaled as Eq. ([7]) and the 
eigenf unctions are distorted. Because the Zeeman field 
h is not scaled, in contrast to the vector potential A, 
the Maki parameter a changes from a as expressed in 

Eq. (I32D. 

We expand the gap equation with respect to the 
operators t), and obtain the eigenfunctions 



(41) 



which are indexed by n and qz. The upper critical field 
equation is decoupled into those for each eigenfunction 

as 



log( 



T . 
r(0)' 



ttT 



dt 



1 







sinh(7rTt) 



r/2 







-1 



1 — cos{ht) cos[—qzVFtcos6'\ exp[ 



m=0 



(rn\Y{n — m)\ 



(42) 

The physical upper critical field Hc2{T) is the highest 
solution of H among the solutions of Eq. at each 
fixed T . In other words, the parameters n and qz are 
optimized so that Hc2 is maximized. 

In systems with A = 0, the pure FFLO state oc- 
curs as described in section IIIII at low temperatures. 
When the system is isotropic, there is infinite degen- 
eracy with respect to the direction of q. We express 
the q's that give the maximum FFLO critical field as 
9 — (go sin0 cos iy9, go sin0 siniy9, go cos 0) with an opti- 
mum value of go and arbitrary and Lp. 

In an anisotropic system with rrix 3> my , vriz , the effec- 
tive Maki parameter becomes large so that a ^ a from 
Eq. ([32]), and S,h > in from Eq. HH). When m^, is very 
large, we can make approximate |<7| « go- Therefore, 
from Eqs. ([57)1 and ([55]) . we obtain the optimum value of 
qz, as 



/ 2n , / 2 2na 



(43) 



for a given n. If we consider that n is always finite for the 
limit a ^ oo, Eq. ([15)) is reduced to gz = ±goj and the 
-ffc2 equation ([42)) is reduced to that of n = 0. Therefore, 
when m^, is very large, states with any finite n have upper 
critical fields very close to that of the state with n — 
and qz — ±go. However, because the index n needs to 
be optimized for each situation, it can be infinitely large. 
For the state with infinitely large n, such that n oc a cx 
(fj, Eq. ([43)) results in jg^j < go and gj_ 7^ 0. 

This can be verified also by numerical calculations. 
Figure [1] shows the temperature dependences of the criti- 
cal fields of states of various n. At each temperature, the 
state with the highest critical field is physical. It is found 
that vortex states with n = 0,1,2 occur depending on the 
temperature, and that the envelope is very close to the 
curve of the two-dimensional limit. It is also found that 
the wave vector q becomes nonzero belowTw 0.51 xTi°^ 

Figure[2]shows the behavior of the critical field for each 
n at low temperatures. We obtain ^ below the tem- 
perature where the solid and dotted curves branch off. 
Interestingly, below T 0.18 x TjP\ the upper critical 
field of the coexistence state with n = 2 and g^ 7^ ex- 
ceeds that in the limit Zm 00. We will discuss this 
later. 
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FIG. 1: Temperature dependences of the upper critical fields 
for Zm = The solid curves show the upper critical fields for 
n = 0, 1,2,3,4,5. At each temperature, the highest field is 
the physical result of the critical field. The broken and dotted 
curves show the upper critical fields in the two-dimensional 
limit and that with the assumption of q = 0. 

Figure [3] shows the Zm (oc a) dependences of the upper 
critical fields i/c2 for «- = 0,l,---8. It is found that 
the upper critical fields of all the coexistence states with 
n 7^ tend to approach that of the coexistence state 
with n = 0, and slightly exceed it, where Zm is large. 
At each 2:„i, the highest upper critical field among those 
with n = 0, 1, 2, • • • is the physical upper critical field. 
It is found that n of the physical state increases as Zm 
increases, and the physical critical field is larger than the 
critical field of the n = state for z^ ^ 1-3. 

Figure ID shows the behavior of Hc2 for large z^. Af- 
ter reaching maxima, all the critial fields for n 7^ de- 
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FIG. 2: Temperature dependences of the upper critical fields 
for Zm = 3 at low temperature. The solid curves show the 
upper critical fields for n — 0,1,2,3. At each temperature, 
the highest field is the physical result of the critical field. The 
dotted curves show the critical field when q = is assumed. 
The thin broken and dotted curves show the upper critical 
field in the two-dimensional limit and that with the assump- 
tion of q = 0. 




FIG. 4: Zm dependences of the upper critical fields for n — 
0, 1, 2, ■ ■ • 8 at T/Tc = 0.01 in the large Zm region. The dotted 
curves show those when q = is assumed. The thick and 
thin curves show the results for n = 0, 1, 2, 3, 4 and those for 
n = 5,6, 7, 8, respectively, a ^ 7.39 x Zm. 



VII. SUMMARY AND DISCUSSION 



crease, and converge with the curve for n = in the 
limit Zm oo. Within the present theory, the physi- 
cal upper critical field, i.e., the highest field at each Zm, 
decreases as Zm increases for very large Zm ^ 7.7, i.e., 
^/zm ^ 0.13. This result seems inconsistent with the 
naive expectation that a reduction of the orbital effect, 
for example, by increasing v-p should cause a weaker pair- 
breaking effect. Presumably, for such a large Zm, a first- 
order phase-transition to a coexistence state must occur 
at a higher critical field than that obtained here, and the 
resultant critical field must monotonically increase with 



a:" 

51 




FIG. 3: Zm dependences of the upper critical fields for n = 
0, 1, 2, ■ ■ ■ 8 at T/Tc = 0.01. The dotted curves show those 
when q = is assumed, a ~ 7.39 x z,n. 



We have examined quasi-two-dimensional type-II su- 
perconductors and the two-dimensional limit. When 
H II [0,0,1], the effective Maki parameter a is propor- 
tional to {rrixmy /miy^^ . Therefore, when nix 3> m^, 
the orbital pair-breaking effect becomes weak, and the 
superconductivity survives up to a higher field, where 
the FFLO state is favored. 

In the coexistence state, the FFLO modulation occurs 
in the direction of the magnetic field, and is smoothly 
reduced to that of the pure FFLO state in the two- 
dimensional limit {mx/mz — + 0), when the directions of 
the magnetic field and q of the pure FFLO state coin- 
cide. In contrast, when their directions differ, it may 
appear that the coexistence state is not reduced to the 
pure FFLO state continuously. However, in actuality, 
modulation perpendicular to the magnetic field is real- 
ized in large n vortex states. The physical origin of 
the order parameter modulation in vortex states with 
higher Landau levels is the spin polarization energy as in 
the pure FFLO state. The coexistence states with opti- 
mum n have upper critical fields close to that of the pure 
FFLO state. Hence, a cascade transition occurs when a 
is large, analogously to the exactl y tw o-dimensional sys- 
tem in a tilted magnetic field ^,m, 25, 26, 27, 28, 29|. 
In the two-dimensional limit rrix —>■ oo, the Landau level 
index n increases as rt cx a oc rrix^^. As a result, the 
coexistence state indexed by (n, gy) (or the pure vortex 
state indexed by n when gn — 0) is continuously reduced 
to the pure FFLO state with q = {q\\ ,q±). 

The relations (|37p and ([55]) connect the coexistence 
states of Eq. (|^T|) and the pure FFLO state in the two- 
dimensional limit. The pure FFLO state with q± ^ 
corresponds to the coexistence state with n on a ^ oo, 
and the upper critical fields of the coexistence states con- 
verges to that of the pure FFLO state. This behavior has 
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been confirmed also by numerical calculations. 

From these results, we can conclude that the FFLO 
state obtained in a theoretical model without orbital ef- 
fects may emerge as the vortex states with higher Landau 



level indexes in real materials where orbital effects are in- 
evitable. In particular, the order parameter modulation 
due to the higher Landau level index is a mark of the 
FFLO modulation perpendicular to the magnetic field. 
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